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Cm Abstract. The perturbative master equation (Bloch-Redfield) is extensively 

» I used to study dissipative quantum mechanics — particularly for qubits — 

C^ despite the 25 year old criticism that it violates positivity (generating negative 

^H probabilities). We take an arbitrary system coupled to an environment containing 

^H many degrees-of-freedom, and cast its perturbative master equation (derived from 

1 ^ a perturbative treatment of Nakajima-Zwanzig or Schoeller-Schon equations) in 

^. 1 the form of a Lindblad master equation. We find that the equation's parameters 

^ are time-dependent. This time-dependence is rarely accounted for, and invalidates 

^_^ Lindblad's dynamical semigroup analysis. We analyze one such Bloch-Redfield 

(-H master equation (for a two-level system coupled to an environment with a short 



Ph 



X 



but non- vanishing memory time), which apparently violates positivity. We show 
analytically that, once the time-dependence of the parameters is accounted for, 
positivity is preserved. 
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1. Introduction 

No system is truly isolated from its environment, thus all quantum systems experience 
some amount of dissipation and decoherence [D H]. To understand the properties 
of real quantum systems we must understand the effect of dissipation in quantum 
mechanics. This is extremely relevant to recent works on qubits and quantum 
information processing (quantum computing and communication). In experiments 
[21 m [5] the coupling to the environment is typically not as small as would be required 
to build a quantum computer. One must understand the effect of the environment on 
a qubit, if one wishes to minimize it. 

Any theory for a quantum system which exchanges energy and information (but 
not particles) with its environment should give a master equation (evolution equation) 
for the system's density-matrix which satisfies three basic requirements; 

(i) preserves the Hermiticity of the density-matrix, so all probabilities are real, 

(ii) preserves the trace of the density-matrix, then the sum of probabilities over any 
complete set of orthogonal states is one, 

(iii) preserves positivity. A system is positive only if the probability of all possible 
states is positive. Given (ii), this guarantees that all probabilities lie between 
zero and one. In this work we do not consider complete positivity, excepting 
comments in Sections |2] and [T] 

There are only a small number of models for which such master equations can be 
derived exactly (we will not address these here). In all other cases, there are two main 
methods for finding such a master equation j2j; 

• Phenomenological method. Here one attempts to construct general master 
equations which satisfy requirements (i-iii). Under the assumption that the 
evolution is translationally invariant in time (a dynamical semigroup property), 
as is often the case for Markovian evolution, Lindblad |7][S] considered the master 
equation given in Eqs. (la|16). He proved that it is the most general equation that 



satisfies (i-ii) above, while also preserving complete positivity. Complete positivity 
is as strong or stronger than positivity, thus it automatically satisfies (iii) above 
(see the comment in Section [t] due to [6]). 

• Perturbative method [1]. Here one takes the evolution of a system and its 
environment (from their combined Hamiltonian) , and traces over the environment 
degrees-of- freedom. Various methods of doing this exist; Bloch-Redfield[51 [TU]. 
Nakaiima-Zwanzig|lll [T^. Schoeller-Sch6n[T3]. However one is typically forced 
to treats the system-environment interaction perturbatively, then all these 
approaches reduce to Bloch-Rcdfield's. 

The Lindblad master equation (the most general generator of a dynamical semigroup) 
takes the form; 



dt 



m= -i[Ksys,m]^- E Y(4^n/5(i)-Hp(t)4^„-2L„p(i)4), (la) 

n=l 

with A„ > for all n, (16) 

where the commutator [A, -B]- = AB — BA, and {iJ„} is a set of ortho- normal (trace- 
class) operators. It is often assumed that all Markovian master equations fall into 
the category of dynamical semigroup evolution, and thus Eqs. ( la|16 ) give the most 
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Figure 1. Timescales for the model (in SectionlGjI fo r w hich we show positivity, 
despite it not satisfying Lindblad's requirement, Eq. ( |16| . The system is a spin- 
half with Hamiltonian —^Ba^, and the environment couples to it via a^. The 
environment's noise spectrum (with noise power So) is broad, leading to a short 
memory time, tm- Decoherence and relaxation times (both ~ Sq ) can be smaller 
or larger than the Larmor precession period, B~^ . 



general Markovian evolution. However this is a subtle point, we discuss it (and define 
terms like "dynamical semigroup" and "trace-class" ) in Section pj 

The perturbative method's advantage over the phenomenological method is that 
one can study how a particular environment (with a given spectrum, temperature, 
etc) affects the system. Thus one can address a crucial aspect of qubit research; 
how should one engineer a particular system to minimize decoherence? However the 
resulting Bloch-Redfield master equation has long been criticized [14] [15] , because it 



can be written in the form in Eq. ( 1 a I but then typically violates Eq. ( 1 & I . In these 



cases it violates Lindblad's condition for complete positivity. Further, there is plenty 



of evidence that it also violates positivity (see Section 1.2 1 



1.1. Outline of this article 

The objective of this article is to study this apparent contradiction between the 
perturbative method and Lindblad's proof. We start by discussing, in Section [2j 
the assumptions that underlie the Lindblad master equation. In Section [3] we consider 
the Bloch-Redfield equation for an arbitrary system, and show that, in general, one 
coupling constant, A2, is negative. However we also show that the parameters of the 
Bloch-Redfield master equation, {A„} and {Ln}, are time- dependent. This means the 
master equation does not generate a dynamical semigroup. Thus Lindblad's proof is 
inapplicable to the Bloch-Redfield equation, and a priori we do not know whether a 
negative A2 will lead to a violation of positivity or not. 

In Section |6] we consider the Bloch-Redfield equation for a particular system 
(a two-level system coupled to an environment with a very broad spectrum of 
excitations). We divide the evolution into two overlapping regimes; short- and long- 
times (sketched in Fig. fll. The time-dependence of the parameters is only relevant 
in the short-time regime {t much less than decoherence/relaxation times). We show 
analytically that the system remains positive in both regimes (i.e. for all t > 0), despite 
the negative coupling constant, A2. 

1.2. The place of this work in the literature 

In traditional derivations of the Bloch-Redfield master equation [T, "5], it is assumed 
that the parameters of the master equation are time-independent. In reality all 
environment-induced terms in the master equation are zero at t = (defined as the 
time at which the system and environment are in a factorized state), before growing 
with t and saturating at i ^ im, where tm is the environment memory time. So the 
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assumption of timc-indcpcndencc is flawed for times of order the memory time, tm- 
This has been discussed in the context of coupled classical oscillators (TB] , over-damped 
Brownian motion [T^, damped quantum oscillators [THl HHJ [23 HH 112, dissipative 
two- level systems [231 [HI (IHl ^S\ and more generally [5SJ [57] . Nearly all these works 
consider dynamics on times of order the t^ as an initial-slip, after which the dynamics 
is given by the time-independent master equation, the justification for this is sketched 
in [Appendix B[ Of most relevance to us are those works which try to show that 
positivity is preserved in this context [23J [HJ [501 [HI [13 HZj • However these works 
provide only plausibility arguments |28j . or numerical studies (they evolved a finite 
number of initial conditions and checked that negative probabilities did not emerge). 
In contrast, for our model, we consider all possible initial conditions and thereby prove 
analytically that positivity is preserved. 

It has been noted that course-graining can ensure positivity |29j . The work 
presented here indicates that the usual assumption of time-independent parameters 
in the master equation only leads to a violation of positivity for t < tm- Thus course- 
graining on such a scale could hide such a violation. It is also common to simplify 
Bloch-Redfield equations by making a rotating- wave approximation [T4t [T5| [30] which 
is also a form of course-graining since it "averages out" fast oscillations. However if 
we treat the time-dependence of the parameters correctly, the Bloch-Redfield equation 
is derived without any approximations which fail on short timescales, so it should 
preserve positivity without any course-graining. 

There has been a lot of interest in a particular class of non-Markovian master 
equations which are positive by construction. They are either constructed by averaging 
Markovian master equations |3T], or by measurement processes |32j- However, while 
these models are extremely interesting, we are not aware of works relating them to 
microscopic models of a typical qubit experiencing dissipation [33] . 

Finally, we mention that some works suggested that the reason for negative 
probabilities was the choice of factorized initial conditions [50l [34] . They argued that 
this initial condition was unphysical, and a more physical initial condition would not 
generate negative probabilities. However, factorized initial conditions correspond to 
any situation in which one makes a projective measurement of the system state at the 
start of the evolution. Thus, while other initial conditions are worthy of study |351I3S] 
(and highly relevant to certain experimental protocols), a factorized initial condition 
is not unphysical, and thus should not be able to generate negative probabilities. In 
this work we restrict ourselves to factorized initial conditions (see Section pi . 

2. The Lindblad master equation 



The Lindblad master equation, Eqs. ( la|l6l, is written in terms of a set of iV^ trace 



class operators, {-£„} (where N is the number of levels of the system). Operators are 
trace-class if they form a complete orthonormal basis in the space of system operators, 
with the scalar-product defined as (L| ■ Lj) = tr[L|Lj], see (37). The basis is complete 
if any system operator can be written as Ogys = J2i ^j tr jLlOsys] . We choose Lq to be 



proportional to the unit matrix. One can see that Eq. (la I preserves the Hermiticity 



and trace of the system's density- matrix (the latter requires cyclic permutations inside 



the trace) . The combination of Eq. ( 1 & ) with Eq. ( 1 a I guarantees positivity. In fact 
it guarantees a stronger condition called complete positivity, which is the requirement 
that all probabilities remain positive even if the system became entangled with a 
second system at i < 0, but then does not interact with it again. For a review see 
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Sections 2.4 and 3 of Rcf. |2,, section VB of Ref. [15^ or the introduction of Ref. [36] , 
In this article we concern ourselves with studying positivity not complete positivity, 
however it has recently been shown ^ that the two are equivalent for the model that 
we study in section [6] 

Lindblad proved that Eqs. ( |la|16 ) give the most general dynamical semigroup 



evolution [71 18]. However to understand if this is applicable to a given system, one 
must ask if that system has the properties of a dynamical semigroup. For this one 
looks at the density-matrix propagator K(i;io), which acts on the density-matrix at 
to to give the density-matrix at time t, so in terms of matrix elements 

P^'J'it) = '^^i'fvijit'^to)Pij{to)- (2) 

This super-operator, K(i;to), isanA^xA^xA^xA^ tensor which acts on the 



N X N density-matrix. Substituting it into Eq. ( 1 a I gives a master equation for 
M^iiji.ij(t). The requirements for Kit; to) to form a dynamical semigroup are given 
in Refs. [2J |H1 [TH], they include (i-ii) above and complete positivity. However 
another crucial requirement is that the propagator must be translationally invariant 
in time, so K(t;to) = IK(t — to) for all t,to > (where the system and environment 
were in a factorized state at time t — 0). Only then docs K.i"j";ij{t2 + ti) = 
'^j^,j,Ki"j";i'j'(t2)M.iij>;ij(ti). Thus a master equation must have time-independent 
parameters to have this semigroup property. If either the system Hamiltonian or the 
environment couplings (coupling constants A„ or operators i„) are time-dependent, 
then K(i;io) is not translationally invariant in time. Thus Lindblad's proof is 
inapplicable for such systems, even if their evolution is Markovian (in the sense that 
d/5(i)/dt is a function only of p{t) not p{t' < t)). So if A„ or i„ are time-dependent 
(as in our perturbative analysis) one cannot a priori state that negative A„ will lead 
to a violation of positivity. 

3. The perturbative (Bloch-Redfield) master equation 

We assume that the system and environment start (at i = 0) in a factorized state 
p(t = 0) (8) Penv This would be the case if the experiment started with a perfect 
projective measurement of the state of the system |38J. The "universe" (system -|- 
environment) then evolves under the Hamiltonian, 

where T and x are system and environment operators, respectively. We treat these 
operators as Hermitian, because we assume they are observables (i.e. charge, magnetic 
dipoles, etc) as is the case in most qubit experiments (and more generally). Without 
lose of generality we can assume T is dimensionless and x has units of energy. 

For a suitable environment one can derive the Bloch-Redfield master equation for 
the evolution of the system's reduced density-matrix, p(t), from the evolution of the 
universe's state (tracing out the environment at time t) . The assumptions necessary to 
derive this master equation are discussed in [Appendix A[ Broadly speaking one needs 
an environment with a broad (almost) continuous spectrum of excitations, then the 



memory kernel of the environment (defined in Eq. (4c I below) decays on a timescale 
tm- Typically the Bloch-Redfield master equation is valid when the memory time, t 
is much less than timescales associated with dissipation (relaxation and decoherence) 



mi 
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which go hkc l/dTxpini) (we set h — 1 throughout this article). The Bloch-Redfield 
master equation can be written as 

^p(t) - -i[J{,y„p(i)]_ -fSp(t) -p(t)Stf + S/5(t)f + fp(f)St, (4a) 

with r being the operator in Eq. dsl and 

S= / drair) exp[-msysT]texp[iiCsysT]. (4&) 

Jo 

UnHke many derivations we do not assume that we can take the upper-bound on this 

integral to oo. The function a^r) is the environment's memory kernel, given by 

a(T) = trcnv [x exp[-i5{cnvT] X exp[iJ{onvT] Pcnv{t)] ■ (4c) 

Since a(T) is typically complex, S is not usually Hermitian (unlike F). We assume 
that Oi{T) is independent of t, then a(— r) = a*{T). This is true if the environment 
is large enough that it is unaffected by the system-environment coupling (during the 
experiment), and the initial environment state obeys [!Kcnv,Pcnv] = 0. The latter is 
the case if the environment is in an eigenstate or a classical mixture of eigenstates 
(such as a thermal state). We assume that air) is a decaying function of r, and define 
the memory time, im, as the timescale of that decay. Then S is i-dependent, because 



t appears in the upper-bound on the integral in Eq. (46) 



Eq. (4a) looks Markovian, in the sense that the rate of change of p{t) depends 
only on the value of p{t) (not the value of p{t') for t' < t). Despite this memory effects 
are present in the memory kernel, a(r). As we see in Appendix A[ if a(r) is finite for 



a given r it means the rate of change of p at time t is affected by p{t — r). This is the 
reason for the time-dependence of S, which is zero at i = 0, and grows to saturate on 
a timescale of order the environment memory time, im- 

BywritingfSp-pStf=i[(fS+Stf),p]^-i[i(fS-Stf),/5]_,where[A,B]± = 



AB ± BA are the anti-commutator/commutator, Eq. (4a) becomes 

^p(t) = -i[:k^y3, m] _ - \ [(rs + 5^f ), p(t)] ^ + sp(i)f + tpm\ (5) 

where we define IK^y, = IKsys - 5i(f 5 - S+f ). Even when S =^ S^, both (f S + S^f ) 
and i(rS — S^F) are Hermitian. The fact that 'K' is Hermitian means that we can 
interprete it as a renormalized system Hamiltonian. 

It is very convenient to define the symmetrized and anti-symmetrized spectral 
function of the noise, S{ijj) and A(uj) such that 

i[a(r) + a{-T)] = Re[a(r)] = j ^^(a;) exp[-iH, (6a) 

/I 
-^A{uj)eM-iuT], (66) 

remembering that we set h = 1 throughout. One can extract the form of S(uj) and 
A{ijj) from environment details (a bath of harmonic oscillators ^39j, a bath of spins 
|40j . etc). For an environment in thermal equilibrium at temperature T |41j . S{lo) 
and A(u)) are related via A{uj) = S'(cj) tanh(a;/2fcBT) |42]. For harmonic oscillators, 
A{uj) (X J {to) and so S{uj) cx J{uj)cot\\{Lo/2k^T), where J(a;) is the spectral-density 
in Ref. |M1. 
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3.1. Dephasing and Lamb shift when a rotating-wave approximation is reasonable 

When the dynamics is dominated by the system Hamiltonian (ofF-diagonal matrix 
elements decay over many Larmor oscillations), then we can make a rotating-wave 
(or secular) approximation |35] of Eq. ([s]). We write p{t) in the eigenbasis JCsys (so 
^sys;ii = EiSij), then we can expect pl°*'{t) — e'(^'~-^j)*/5sys;ij(i) to be insensitive 
to all fast oscillating contributions to its dynamics. We neglect ("average out") 
contributions to {dp\°^ /dt) which come from Psys;i'j' when i' ^ i or j' ^ j, since 
these contributions oscillate fast, at a rate {Eii — Eji — Ei + Ej) [ll]. Then 
(d/dt)/5™* = [iA_E(i, j) — T2"^(i, j)]p™*. The dephasing rate, at which a super-position 
of states i and j decays to a classical mixture (I/T2 for two-level systems) is 

T^\i,j) 0, Re[i(fS + Stf),, + i(fS + Stf)„- - S,,f„- - Tu^,] • (7) 

The coupling to the environment also causes a Lamb shift; the precession rate is 
modified by the sum of the modification in !K' and AE{i,j), where AE{i,j) is the 
imaginary part of the square brackets in Eq. (l7| . 

3.2. Writing the Bloch-Redfield equation as a Lindblad equation 

To cast Eq. (pi in the Lindblad form, we rewrite it in terms of a set of orthonormal 
(trace-class) operators, {Pi}. We use the usual Gram-Schmidt procedure; defining 
Pi oc r, and P2 as proportional to the component of S which is orthogonal to F. The 
constants of proportionality are such that both Pi and P2 are normalized. Hence 

f 
Pi = 



tr[f2] 

tr[EtHl-|tr[P,'E]|' 

SO Pi is Hermitian while in general P2 is not. As Pi,P2 form an orthonormal basis, 
we have f = A tr[Pir] + P2 tr[P2f ] and S = A tr[As] + A tr[As]. Then the Bloch- 
Redfield equation becomes 

^p(t) = - i[M^y„ p{t)] - \ E h^AHp^m + mplp^ - ^p^mpl)- (9) 

In general, hij = tr[f ^A] tr[A^S] + tr[StA] tr[A^f]. However here tr[Pjf] = 0, so 
hij is given by the ijth element of the matrix 

where for the compactness of what follows we have defined 

6+ = tr[ftA]tr[A^S], 

5, = Re(tr[ftA]tr[As]), (11) 

we also define b^ = |6_|-p-|-6^. We retain fs on the symbols to make the structure clear, 
however F^^ = F and A — A oc F. The eigenvalues, Ai_2, and the SU(2) rotation, U, 
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to the eigenbasis of h, are 

Al,2 = 6^ ± &, 



U = 



1 



(l + 6,/6)i/2 



\A(ft+M 






(12) 



Performing this rotation on Eq. S , the Bloch-Redfield equation takes the form of the 
Lindblad equation, Eq. {la), with Li = J2j=i 2^y-^i 05]- However in general A2 is 



negative [T31[Tn], so this does not satisfy Lindblad's requirement in Eq. (lb) 



4. Perturbative master equation for an extremely short memory time 

We assume here that the memory time, im, is much shorter than any timescale in !Ksys, 
i.e. im ^ ^^i where Agyg is the largest energy difference in the system's spectrum. 
We substitute T{-t) = f-i[ic,y„f]_T- i[a<,y„ [:H:sys,f]_]_T2 + 0[(AsysT)3] into S. 
We expect that a(r) is always given by a dimensionless function of r/i,„ multiplied 
by i~^ (given that h—1). Then S(t) (having units of energy) is 



s(i) = /o(t)r -it„ji(i)[:K,ys,r]_ 



2^mJ^ vy L-"-sysj 



0[t-XA,ydm)'], 



(13) 



where fq{t) — J^ dr {r/t^Y a{T). For all g, fq{t) goes like i~^ multiplied by a 
dimensionless function of i/im- Thus Eq. (13 1 is an expansion to second-order in 
Writing fq{t) in terms of S{uj) and A{ijj) we have 



powers of A 



sys^m- 



IS) 



where a positive infinitesimal constant, 



-(5(.) + A(.)) — 



1 



Aujt 



iO^ 



ensures the convergence for t 



(14) 
00. Thus 
■ (15) 



S(i) - (/o(t)Vtr[r2] - \t1h{t)K)P^ + t„Ji(t)^2tr [rj{sys[?i:sys,r]_] P2. 

The only /2(t)-term that we keep is in the prefactor on Pi, for compactness we define 



K = tT[T['Ksys, [IKsys, r]_]_]/y tr[r2]. This term gives a ©[i^J-term in the final result, 
while other such /2(t)-terms give at worst a 0[i^]-term. From Eq. (11), we get 



b+ = t^ hit) J2tr [r2] tr [rj{3y.[J{sys,r]_] , 



b, = Re[/o(<)] tr [r^] - ^^tlRc[h{t)]KJtT [T^] 



(16a) 
(166) 



Since tr [fjCsyslIKsys.r]-] < A^^^tr [f^], we have \b+\ - {A^y^t^)b^ where A^ysi^ < 
1. The terms that we dropped only give contributions of order (Asysim)^^^. 

To zeroth order in im we recover Lindblad's result, Eqs. ( la|16 l, with only one 
non-zero coupling constant Ai = 2Re[/o]tr [F^] > 0, associated with the operator, 
Li — (tr[r^])~^/^r. However to first order in i,n, we have Eq. [la) with two non-zero 
coupling constant Ai and A2; the latter of which is negative (even for infinitesimal tm)- 
We use this model to explore the contradiction between Bloch-Redfield and Lindblad. 
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4-.1. Environment with a nearly white-noise spectrum 

Here we consider an environment with a nearly white-noise spectrum of excitations (a 
very wide Lorentzian), at extremely high temperature, k^T ^ uj^, so 



S{lli) = 5*0 



W' 



Aico) 



2kBT < 



(17) 



where A{u!) is given by the result below Eq. {6b\. The high-energy cut-off, Wm, plays 
the role of the inverse memory time, i~ ^ , so for nearly white-noise we need it to be 
much larger than the largest system energy scale, Asyg. Then Eq. (14 1 gives 

/o(i) = ^So (1 + i(2fcBTinO"') (1 - cxp[-t/t„,]), (18a) 

/i(i) = iS'o (l + i(2fcBTt„,)-i) (l-(l + t/i„,)expHA„,]), (18&) 

/2(i) = 5o {l + i{2kBTt^)-') (l-(l + tA„, + (tA„,)V2)expH/im]), (18c) 

where we evaluated the w-integrals using complex analysis (by pushing the contours 
into the upper-half plane, one finds that the results are due to the pole at w = iwm)- 
Both fo{t) and /i(i) go exponentially to their long-time limit (i ^ tm), with the rate 
given by the memory time, im- When t/im 3> 1 we have fi(t)/ fo(t) ~ 1, while when 
t/tm <C 1 we have fi{t)/fo{t) — t/t^. For such an environment, the Bloch-Redfield 
equation is valid for S'otm <C 1, see [Appendix A) 



5. Positivity as a constraint on a two-level system's purity 

To ensure that there is no basis in which the density matrix has negative probabilities 
(i.e. no possible measurement will return an unphysical probability) it is sufficient and 
necessary that the density matrix's eigenvalues, {Ak}, satisfy < A^ < 1 for all k. To 
see this, consider an arbitrary basis which is related to the eigenbasis by the unitary 
transformation, 11. In this basis all probabilities are given by pa — ^j, lUi^pAj., where 
the unitarity of 11 guarantees that J2k \'^ik\'^ = 1- Thus if < A^ < 1 for all k, then 
probabilities in this arbitrary basis, satisfy < p^i < 1 for all i. 

A two-level system is special because the eigenvalues of its density matrix are 
defined by a single parameter, s (remember that the sum of the eigenvalues must 
be one). The most general two-by-two density-matrix is of the form p = |(o'o + 
Sx^x + Syffy + SziTz), where ^x,y,z are the Pauli matrices, and Sx,y,z are real numbers, 
when diagonalized it takes the form pd = i('5'o + scTz) with the single parameter 
s^ — s^ + Sy -t- s^. Thus to ensure positivity we require that — 1 < s ^ 1- The purity of 



pisP = tr[p2]^i(l + s2 



j), thus ensuring positivity is equivalent to ensuring 



that P < 1. This is not the case for systems with more than two levels [46] 



Finally it is worth noting that Eq. (|la| leads to 



di 



= 2tr 



m 



dm 

dt 



= -2 



E 

n=0 



A„ tr 



Li[u,,p{t)\^m]- (19) 



6. Two-level system with nearly white-noise: proving positivity 



We now consider a two-level system with "K^ 



— ■^B&z, coupled to an environment 
viaF = (Ta;. Then Eq. (4&I gives S = J^ dTa(T) [a r^ cos Br — ay sin Br] , so Pi = ax/V^ 
and P2 — <^y/V2- For an environment with a short memory time, Eqs. (16a 16&I give 
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Figure 2. A sketch of the Bloch sphere, for the situation discussed in SectionlH] 
In (a) we show the axes associated with j^sys, Pi,2 and Li^2- We show -£1^2 
for a case where they are Hermitian {w is real), as only Hermitian operators are 
associated with axes in the Bloch sphere. In (b) we sketch the effect of the Li,2- 
terms on the evolution of the Bloch vector which represents the density-matrix, 
i"Bloch = (2Re[pi2], — 2Im[pi2], pii — P22)- The Li-term reduces the magnitude of 
the vector in the plane perpendicular to Li (diagonal cross-hatching) at a rate 
given by Ai. The L2-term increases the magnitude of the vector in the plane 
perpendicular to L2 (vertical cross-hatching) at a rate given by IA2I <C |Ai|. 



b+ = 2Bt^ /i(t) and b, = 2Re[/o(t)] - {Bt,nfRe[f2{t)]. Thus to second-order in Bt,„, 
Eq. ^ gives 

,2 f l/iWP 



Ai = 4Re[/o(t)] + (Bt„ 
{Bt^r\Mt)\' 



VRc[/o(i)] 



2Re[/2(t)] 



Re[/o(t)] ■ 
Defining w = i?t,„/i(i)/Re[/o(t)], the Lindblad operators, -£1,2, are given by 



L2 



= U 



Pi 



i-lkP 



w 



Ikp 



(20) 



(21) 



Here we give U to first order in Bt^, but keep the higher order terms necessary to 
ensure U^li — 1. Note that Li and L2 are not Hermitian unless /i(i) is real. 

We take the T ^ 00 limit of the nearly white-noise in Section [4A] so Im[/i(i)] = 0, 
then w is real and gives the angle marked in Fig. [2^. Defining the x', y'-axes such that 
-£1 = o'x' /V^ and L2 = (Ty' /^/2, the Bloch- Redfield equation reduces to 

'^ m = - i[:H;^y„ p{t)] - 2Ai [m - ^.-m^^x') - 2A2 (^m - ^ymc^y) ■ (22) 



di 



The coupling constants, Ai, A2, are given by Eqs. ( 18a|186(20 ) with T — > 00, so 
Ai = 2S'o(l-e-*/*») + 0[(5<m)'], 



A2 — — (Stm) S( 



2, (l-O + VtmK!^ 

2(1 -e-*/*™) 



0" 



(23) 



Substituting these results into Eq. ( [19| , and writing — A2 as +IA2I to emphasis that it 
tends to increase the purity, we get 

dP 

~dt 



-Aitr [p^t) - (a^'mY] + |A2|tr [p^{t) - {ay,p{t)y 



(24) 
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For times, t, much less than 5*0 we can get the purity to first order in 5*0, by 
integrating Eq. ( 24 1 with p{t) replaced by its value to zeroth order in 5*0 



p^'^'{t) = 5 1 + {sx cos Bt + Sy sin Bt)ax> + {sy cos Bt — s^ sin Bt)ay' + s^Oz 



(25) 



where the constants (sa;, Sy,Sz) define p{t — 0). Note that we have used the fact that 
to zeroth order in 5*0 we have Or{'^^^ — 'Ksys — —\Baz. Then Eq. (24 1 becomes 

dP(t) 



di 



= -\i[{syCOsBt~ s^sinBtf + si] + \X2\[{sx cos Bt + SySinBt^ + sf\ 



+0[Slt]. (26) 

As ini ^ B^^, we can restrict ourselves to times t ^ B^^^S(^^ (and hence expand 
in powers of Bt and S^t), and still study the dynamics up to times 3> im- The 
problematic coupling constant, A2, is 0[i?^], so we must expand the right-hand-side 
of Eq. (26 1 to 0[-B^], to see the effect of A2 on the dynamics. After this expansion in 
B, we expand the purity about P(0) = 1. So P{t) = 1 + /J di'(dP(t')/di') gives 



P{t) = 1 - 25otm [si h +slly- 2SxSy Bt,, hy + sl{Btn,fL 

where Iz, ly, Ixy and Ix are the following functions of t/t-^, 
■ ,,Ai(t') + A2(0 



0[St,t,SQBH'] (27) 



/.(t/im) = 

Ix{t/t„,) = 



2oolm 

'^^, a~iBtr)Xi{t') 

2Dotm 

' {BtrXi{t') + \2{t') 

2S0BHI 

t/tyn 



t/t„ 



t/tn 



t/t„ 



dj/(l-e~''), 



dj/(l 







^^(1-e-^)- 



di/i/(l 



[1 - (1 + ^)e-'^]- 
4(1 -e-"^) 



(28a) 
(28&) 
(28 c) 

(28rf) 



where v = t' /t^. The "~" indicates that we keep only the leading order in (Bt^) in 
each term, this will be sufficient for our purposes. 



To show that P{t) does not exceed one (in the range of t for which Eq. ( 27 1 is 
valid), we show that the square-bracket in Eq. (pTJ) is never negative. Writing the 
square-bracket as [izsl + Iy{sy - Sx{Btni)Ixy/Iyr + sl{Bt,^)^{Ix - lly/Iy)], we see 
that there are three terms; the first two are always positive (but will be small for spins 
starting close to the z'-axis, i.e. Sy, s^ <C 1), the third term is positive if Ix > Ixy/^y 
Thus we must show that lylx/Ily > 1- For t <C tm, 



W/traf 



Ixy * 3('-/*mj 



(iA,n)' 



and for t :^ t„ 



I,, 



Thus for t <^ ^m wc have lyl^/I^y - 



W/trnf 



W/tm)' 



(29) 
(30) 



135/128, while for t ^ ini we have lylx/Ixy 



4/3. For finite t we see that lylx/I^y is a monotonic function of t which goes from 
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r(o°) f(°°) //'r(oo)\2 



Figure 3. Plots oi lylx/I^y (solid curve) and ly 7i /{Ixy ) (dashed curve) 
as functions of t/tui- The two horizontal lines are the tvifo extrema of lylx/lxy' 
its small t limit of 135/128 and its large t limit of 4/3. The crucial point is that 



(oo) ,((X)) ,/t(oo)-^2 



which 



h^x/lly > 1 for all t > 0. This is not the case for li, 'li ' /{liy') 

one would get if one mistakenly assumed time-independent coupling constants 

Ai(oo), A2(oo); this is less than one for all t/tm < \/3 and goes to — oo at i = 0. 



135/128 to 4/3 (see Fig. [3]), thus it is always greater than one. This means that 
P{t) < 1 for all times much less than Sq^ (including times greater than tm). 

For t„i < t < B-\5'o"\ the purity P{t) = 1 - 2SQt[sl + {sy - ^! 



0[BH 



3+31 



j^^isxBty]. By checking all pure initial system states (state with s^ 
we see that this is maximal for Sx = i[i + {Bt)'^/8] 
(where Sx and Sy have the same sign) [47]. Hence 

P(t„,«t«B-i,^o-i) < 1 



,Bt} 

^2 _ 



1) 



±^Bt and s^ = 



1 c r2 .3 



(31) 



This upper-bound on the purity will be crucial for our proof (in section 6.2) that the 
purity does not exceed one at longer times. 

// we had made the usual assum,ption that one can replace \i{t) and X2{t) with 
\i{t = oo) and \2{t = oo) for all t, then we would have Eqs. ( 28a|28d l with ly, I 

(oo) _ rt/*m .„ Aoo) _ ft/t^^,,,, _^ r(oo) _^Jtji 



xy 7 ^x 



replaced by Iy°°' = /q " d:/, r^' = J^^' "" dvv, and Ix = /q''" di^ \u 



1/4]. We 

plot ly Ix / {Ixy Y in Fig.bl and see that it goes to — oo as i/tm -^ 0. Thus such a 
mistaken assumption would have led us to conclude (as other have) that P can become 
bigger than one (at times < t,n). The mistake is most clearly illustrated by looking at 

Eq. (27 1 with Sx — ^ and Sy — Sz — 0, then using Ix in place of Ix would lead one 
to think that for i ^ i,„, the purity would be 1 -I- \Sot when the correct expression 
shows it is 1 — |S'otm(i/^m)'^- Thus it is only by keeping the time-dependence of the 
coupling constants, that we can show that the purity cannot exceed one for all times 
<C Sq^ (including times greater than t^). 



6.2. Positivity at long times (times of order and greater than Sq ) 

We now turn to the evolution of the purity at all times much greater than t^ (the 
long-time regime in Fig. 111). For times of order and greater than Sq we need the 
full Bloch-Redfield equation, Eq. (22), not just the short time expansion of it. Since 
t ^ tin, the coupling constants have saturated at their long time limits; Ai = 25*0 and 
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A2 = -2So{Bt^/2f . Then Eq. Q reduces to 

^ = -25o[(l - {Btj2f)sl{t) + 4(t) - {Btj2fsl{t)]. (32) 

Since Bim *^ !> we can see that P{t ^ im) decays for nearly all Sx^y.z{t)- However 
the purity may grow if Sy{t) ~ Sz(i) ~ 0[(Sim)^]; then the purity might exceed one 
(particularly if s.^ is close to one) . 

To see if the purity can exceed one, we expand the evolution about the time t^, 
where we choose io such that Syito) = 0. We then perform the same expansion about 
t — Iq here as we performed about i = in Section |6.1[ Hence on the right-hand-side 



of Eq. (32) we make the substitution Sj:{tQ + t) = s'^ cos Bt, Sy{tQ + t) = s'^ sin Bt, 
Sz{to + t) = s'^, where we define s'^ = Sx{to) and s'^ = Sz(io) (remember that to is 
chosen such that Sj,(to) = 0). This substitution is good for all r ^ Sq^ , After the 



substitution we expand the right-hand-side of Eq. (32 1 up to second order in Bt. Thus 
ior T <$: B-\ S^\ 

^ = ~2So [(1 - (i?t„./2)^)4' + {s.Brf - (45^/2)^] . (33) 

"'- t=ta+T 

From this we see that the purity can only increase during a time-window where 
|t| < [(im/2)^ - (4/sL-B)^]^/^ (neglecting a term that is higher order in Bt^). The 
maximum possible time for this growth is t„^ (i.e. when s'^ ~ 0, P grows during the 
time-window from r = — |im to r = |im)- Thus the assumption that r <C B^^, Sq 
is fulfilled for all r at which the purity is growing. 

At this point it is sufficient to make a gross over-estimate of the amount by which 
the purity can grow. If we assumed that the purity grows during the entire time- 
window — 2^m ^ T < 2^111 at the maximal possible rate (i.e. the rate at r = when 
s^ = and s^. = 1), then during this time-window it would grow by ^SoB^t"^. If we 
define AP as the true increase of the purity in the time- window where it grows, the 
over-estimate enables us to put the following upper-bound; 

AP < ISoB^l. (34) 



Comparing this with the upper-bound on the purity in Eq. (31) with t 3> im (but 
t <$:. Sq^), we see that increasing the purity by AP cannot cause it to exceed one. The 
short- and long-time regimes overlap (see Fig. [I]) , so by showing that P < 1 in both 
regimes we have shown positivity for all i > 0. 

7. Conclusions 

The Bloch-Redfield master equation for an arbitrary system can be written in the 



form of a Lindblad master equation, Eq. (la). Only by setting the memory time equal 



to zero (strictly Markovian evolution) do we recover Lindblad's result with coupling 



constants, {A„}, which are time-independent and positive, Eq. (lb\ 



For finite memory times, the Bloch-Redfield master equation can still be cast in 



the form of Eq. ( |la| , but its do not satisfy Eq. (1&). However, the parameters are 
time-dependent which means that the semigroup property is absent, and so Lindblad's 
requirements are inapplicable. We show analytically for a particular model (a two-level 
system coupled to a high-temperature environment with a memory time much less 
than system timescales) that the master equation preserves positivity if and only if we 
keep the time-dependence of the parameters. Further, it turns out that positivity and 
complete positivity are equivalent for this particular model [B] . 
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It is remarkable that our result only coincides with Lindblad's for strictly zero 
memory time, t^^ — 0. If we take the limit t„^ -^ 0, we find that one coupling 
constant tends to zero from below. Further, we argue (see the appendix) that the 
Bloch-Redfield equations become exact in this limit. Thus even for infinitesimal <m, 
one coupling constant is negative. Positivity (and hence complete positivity) is none- 
the-less preserved by the time-dependence of the coupling constants at times of order 
the infinitesimal time t^. 

We wonder if an analysis of the time-dependent parameters of an arbitrary Bloch- 
Redfield master equation would show positivity, or even complete positivity. If this 
could be proven, one could argue that the Bloch-Redfield master equation contains 
both the Lindblad equation and finite memory-time corrections to it. 
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Appendix A. Deriving Bloch-Redfield from a Dyson equation 

For completeness, we sketch the derivation of the Bloch-Redfield master equation 
(HI EH] , using a common "modern" approach [IT] based on a real-time Dyson equation 
[T5] . The derivation is none-the-less equivalent to Refs. [3 HO]- At i = the system 
and environment are in a factorized state (e.g. a perfect projective measurement is 
made on the system at i = 0). The propagator of the system's reduced density 
matrix is K,-jv.,y(i;0) = trenv [(i'|e-'^"""*|i)/5o„v(j|e'*"»"*|j')] , with Eq. ^ giving the 
system's reduced density-matrix at time t. The Dyson equation for K(t; to) (treating 
the system-environment interaction as a perturbation, which we keep to all orders) is 

K{t; 0) = K'y%t; 0) + f dt^ f ' dtiK'^'it; ^2)^(^2; ii)K^''^(ti; 0), (A.l) 

Jo Jo 

where K{t; t') is the propagator including all interactions; K'^y'^(<; t') is the bare system 

propagator (propagating it only under the Hamiltonian IHgys). Since there are no 

interactions after i = in the first term and after time ^2 in the second term above, 

we can trace out the environment at these times. Finally 5](f2; ^1) an irreducible block 

of the propagator (with the same tensor structure as K(t; 0)), it is the smallest block 

for which the system has interacted with o ne or more environment excitations. 



Taking the time-derivative of Eq. (A.l I, and noting that (d/di) /J dt2K''^'' {t; t2)F(i2) 



¥{t) + /„ dt2E"y'(i)K'^y"(i; t2)F(i2) for any F(i2), we get the master equation 

IK(i;0)== -iE"y"(f)K(t;0)+ / diiS(i; ii)K(ti; 0). (A.2) 



dt JO 

?sys 



We have defined E^?",.,, - {i'\:Hsys\i){msys\f), then (d/dt)K^y^(i; 0) = 



-iE^y"(t)K"y^(i;0). To clearly see the non-Markovian nature of Eq. (A.2 1 we can 
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Figure Al. (a) Real time Dyson equation for an arbitrary system. The pair 
of lines with cross-hatching between them are the full propagator, K(f";i'); the 
lines without cross-hatching are the bare system propagator, W'^''(t"; t')\ and the 
lines with solid colour between them are the irreducible block, S(t",t'). Internal 
indices are summed over and internal times are integrated over as in Eq. ( |A.1| |. 
This drawing of the propagators emphasizes that only K^'f,'' .,, -, ,,{t";t') can be 
written in the form A^ll.^l X Bjh.ji. (b) Some lower-order diagrams for the 
irreducible block, S(i",t'), in all cases we integrate u), u)\ etc, over the spectrum 
of excitations. The second-order diagrams are labelled (1) to (4). 

substitute it into Eq. ^ which gives (d/di)/5(i) = -i[:ksys(t), p(t)]+/o diiS(i; ti)p{ti). 
This master equation is exact, our only assumption was that the system and 
environment were in a factorized state at time t — Q. It is formally equivalent to 
the Nakajima-Zwanzig equation [TTl [12]. However it is of little practical use (giving 
no great advantage over standard perturbation theory) unless the irreducible block, 
S(i2;ii) is reasonably local in time, i.e. decays on a scale t2 — ti <C t. Without 
approximation we can use p(t) = K(t; ti)p(ti) to write this master equation as 

-p(t)= -i[J{,y,(t),p(i)]+ / dtiS(t;ii)K-i(t;ii)p(t). (A.3) 



This might "look" Markovian, but the non-Markovian nature is in the new term 



^{t;ti). Approximations of Eq. (A.3 1 will give a Bloch-Rcdficld master equation. 



Appendix A.l. The Block- Redfield equation from a Born approximation 

Here we get the Bloch-Redfield master equation by making a Born approximation of 



the irreducible block S(t",t') in Eq. (A.3 1. It involves neglecting all contributions to 



S(t",i') beyond second-order. Our derivation involves two assumptions which justify 
the Born approximation (other derivations may be possible). 

Our first assumption is that the environment is large enough to have a continuous 
energy-spectrum of excitations (although it does not matter if this spectrum is discrete 
on scales -^ t^^). So for finite relaxation/decoherence rates, we assume the coupling 
to each environment excitation is small enough to be treated only up to second order. 
Thus each excitation evolves only under Mcnv up to the time of its (first or second 
order) interaction with the system. It then never interacts with the system again, so 
we trace it out immediately after the (first or second order) interaction. 

Our second assumption is that the environment's initial density-matrix obeys 
[IKcnv, Ponv] — 0, as would be the case for either an eigenstate or any classical mixture of 
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eigenstates of Jfonv (such as a thermal state) . Combining this with our first assumption 
means that we can treat ponv as time-independent. Then without loss of generality 
we can make tronv [ipcnv] = 0, by moving any constant off-set into the definition of 
JCsys- This removes the first order contributions from the irreducible block, 'E{t",t'). 
Thus S(t",t') becomes the sum of second-order (and higher-order) terms sketched in 
Fig. |Al[ 3. The dotted lines indicate that a given environment excitation (with energy 
iv) has been created by the system-environment interaction. 



Treating the integral in Eq. (A. 3 1 to lowest (second) order in x, means making 
a Born approximation on S(t",i') and treating K^^(t;fi) to zeroth order in x [45]. 
Hence defining r = t - ii and S(t) = SB°™(t; t - T)[Wy%t; i - r)]-\ we have 



^p(t) = - i[^,y,{t),p{t)] + J drS(T)p(i). (A.4) 



The four contributions to S^°™(t; i - t), labelled (1-4) in Fig. \Alp , give 

4'i';.,(^) = tW [(*'|f xe-'^«-f xe'*''n*>Ponv(j|j')] , (l-5«) 

4?';„-(^) = tw [(^'|^)penv(J|e-'*°^f ie'*°^f x|j')] , (1.56) 

^k^J(^) = ti-^v [(*'|e-'^''^f xe'*«n*)pe„v(j|f x|/)] , (1.5c) 

4?';..(^) = t^^v [(z'|f x|*)penv(j|e'*"^f xe-'^«n/>] , (l-5d) 

where we define JCq = Jfgys + 3^^onv We re- write all these contributions in terms of 
operators acting to the left and right of the density-matrix, p{t). Those interaction 
on the upper line are to the left of p{t), while those on the lower line are to the right. 
Thus summing these four terms we get S{T)p{t) = tronv [r(0)a:(0)r(— t)x(— T)[p(t) (g) 

Pc„v]+[/5(i)0Pc„v]f(-T)x(-T)f.T + f(-r)f(-T)[p(t)0pe„v]f(O)i(O)+f(O)i(O)[p(t)® 

penv]r(— T)a;(— r)] , where the operators are in the interaction picture, so T{t) = 

exp[iIKsysT]rexp[— i!KsysT] and x(t) = exp[iJ{env''']2;exp[— iJCenv''"]- Substituting this 
into Eq. (A.4 1 we get the Bloch-Redfield master equation that we gave in Section [s] 



Finally, to see when the Born approximation is justified, we must estimate the 
higher-order contributions that we are neglecting. The higher order contributions to 
S(i; ii) take a similar form to the second-order ones, but have more factors of Fa; acting 
to the left and right of the density-matrix. The times at which these interactions can 
occur are chosen such that that E(i;ti) is irreducible (as d iscu ssed above). A typical 



fourth order contribution (those in the second line of Fig. Al d) will go like |Fa;| i 



14+3 



^ini 



compared with the second-order terms that went like iFaJptm. It is justifiable to neglect 
the fourth-order while keeping the second-order, only if |ra;|ti„ <C 1. Physically the 
constraint that Irxjim ^ 1 means that the Bloch-Redfield master equation applies to 
situations where the decay rate of memory effects, i~^ , is much faster than dissipative 
(relaxation and decoherence) rates ^ iFxpim- There is no constraint on the ratio of 
dissipative rates to the system's energy-scales, so the Bloch-Redfield equation can be 
applicable to strong (over-damped) and weak (under-damped) dissipation. 

Appendix B. A simple picture of initial-slips 

To understand how initial-slips work [Tni[T71[IHlllSl[lSl[inil21[ini[2IJIlS],it helpful 
to neglect the matrix structure of the master equation. Then one has 

{d/dt)p{t)^~F{t)p{t), (2.1) 
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where F{t) is time-dependent, but saturates at a finite value, fac, for times greater 
than the memory time {F and p are now numbers not matrices). This is traditionally 
approximated by [TJ |5] 

(d/di)p(i) = -/ooP(t). (2.2) 



Eq. (2.2) gives the wrong evolution for any initial condition, p(0), because it is not 



justified for times less than the memory time. However, by multiplying p(0) by an 



initial-slip one can ensure the evolution under the incorrect Eq. (2.2) coincides with 



the evolution under the correct Eq. (2.1 1 for all times much greater than the memory 
time. For the above equations, the initial-slip is simply exp [ — /q dt'[F{t') — /oo]] . For 
t much greater than the memory time, tm, the initial-slip becomes time-independent 
(one can take the integral's upper-limit to oo since [F{t') — /oo] ^ for t' :S> im)- Thus 
one can take p{0), "slip it" so that it becomes exp [— J^ dt'[F{t') ~ /oo]] /o(0), and use 
that as the initial condition for evolution under the incorrect Eq. ( |2.2[ ) . The resulting 
p{t) coincides with the correct result for all times much greater than t^, but will 
be absolutely meaningless for all times of order t^. Qualitatively the same analysis 
applies to the full master equation, but it is complicated by the matrix structure of 
the master equation (see e.g. Ref. jM]). 

The above sketch of the initial-slip method, makes it clear that it is not suited 
to our analysis of positivity. The short-time dynamics (on timescales of order the 
memory time) that it generates are fictitious; a sudden initial-slip of the density matrix 
followed by evolution under an incorrect master equation. Studying the positivity for 
these fictitious short-time dynamics tells us nothing about whether the true short-time 
dynamics preserves positivity or not. 
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